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Abstract 

We investigate the differential branching ratio, branching ratio, differential forward-backward 
asymmetry, the forward-backward asymmetry of the lepton pair and the lepton polarization asym- 
metry of the exclusive B —>■ t^t^^ decay in the general two Higgs doublet model including 
the neutral Higgs boson effects. We analyse the dependencies of these quantities on the model 
parameters and also on the neutral Higgs boson effects. We found that they get considerable en- 
hancement from the two Higgs doublet model compared to the standard model and neutral Higgs 
boson effects are quite sizable. 
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1 Introduction 



Rare B-meson decays are induced by the flavor-changing neutral currents (FCNC) and they occur 
only through electroweak loops in the standard model (SM). Thus, on the one hand, they provide 
fertile testing ground for the SM and on the other hand, they offer a complementary strategy for 
constraining new physics beyond the SM, such as the two Higgs doublet model (2HDM), minimal 
supersymmetric extension of the SM [p, ect. From the experimental point of view, studying rare B 
meson decays can provide essential information about the poorly known parameters of the SM, like 
the elements of the Cabibbo-Kobayashi-Maskawa (CKM) matrix, the leptonic decay constants etc. 

In this work, we study the radiative B — > t^t^'J decay in the general two-Higgs doublet model 
(2HDM). It is induced by the pure-leptonic decay B t^t^ and in principle, the latter can be used 
to determine the decay constant fs [@], as well as the fundamental parameters of the SM. However, 
it is well known that processes B £^£^ are helicity suppressed for light lepton modes, having 
branching ratios of the order of 10^^ for i = ji and 10^^° for £ = e channels [j^]. Although the 
B r+r^ channel is free from this suppression, its observation is difficult due to low efficiency. If a 
photon line is attached to any of the charged lines (see Fig.l), the pure leptonic processes B 
change into the corresponding radiative ones, B £+£"7, so helicity suppression is overcome and 
larger branching ratios are expected. Depending on whether the photon is released from the initial 
quark or final lepton lines, there exist two different types of contributions, namely the so-called "the 
structure dependent" (SD) and the "internal Bremsstrahlung" (IB) respectively, while contributions 
coming from the release of the free photon from any charged internal line will be suppressed by a 
factor of ml/M^r. The SD contribution is governed by the vector and axial vector form factors and 
it is free from the helicity suppression. Therefore, it could enhance the decay rates of the radiative 
processes B r^r^'j in comparison to the decay rates of the pure leptonic ones B t^t^ . As for 
the IB part of the contribution, it is proportional to the ratio nii^/niB and therefore it is still helicity 
suppressed for the light charged lepton modes while it enhances the amplitude considerably for ^ = t 
mode. Indeed, B £+£"7 decay have been investigated in the framework of the SM for light and 
heavy lepton modes as well as in the models beyond the SM [7p, and it was found that in the 

SM BR{B -> i+i--f) = 2.35 x IQ-^H, 1.90 x lO^^H, 9.54 x IQ-^ [§, for £ = e, fx, r, respectively. 
With long distance contributions, BR{B r+r"7) = 1.52 x 10"*^ was obtained In 2HDM, in 
contrast to the channels with light leptons, the channel B t'^t^'j receives additional contributions 
from the neutral Higgs boson (NHB) exchanges, in addition to SD and IB ones. In [|^], B t^t'j 
decay is investigated in model I and II types of the 2HDM including NHB effects and shown that 
these effects are sizable when tan /3 is large. Our aim in this work is to study the sensitivity of the 
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physically measurable quantities, such as branching ratio, photon energy density, forward-backward 
asymmetry of the final lepton and lepton polarization asymmetry, to the NHB effects, as well as to 
the model III parameters, like the Yukawa couplings ^^^^ and Cw.rr- 

The work is organized as follows. In section 2, after a brief summary about the main points of the 
general 2HDM, we first present the leading order (LO) QCD corrected effective Hamiltonian for the 
process b st^t^ , including the NHB effects and then give the matrix element for the exclusive 
B T^T^'y decay, together with the explicit expressions for the double differential decay width, 
photon energy distribution, forward-backward asymmetry and the polarization asymmetry of the final 
lepton . Section 3 is devoted to the numerical analysis of the dependencies of these observables 
on the model III parameters and also on NHB effects. Finally, in the Appendix, we give the explicit 
forms of the operators appearing in the Hamiltonian and the corresponding Wilson coefficients. 

2 The B T^T~^ decay in the framework of the general 2HDM 

The 2HDM is the minimal extension of the SM, which consists of adding a second doublet to the 
Higgs sector. In this model, there are one charged Higgs scalar, two neutral Higgs scalars and one 
neutral Higgs pseudoscalar. The general Yukawa Lagrangian, which is responsible for the interactions 
of quarks with gauge bosons, can be written as 

where i,j are family indices of quarks , L and R denote chiral projections L(R) = 1/2(1 =i= 75), 
for m = 1,2, are the two scalar doublets, are quark doublets, Ujr, Djr are the corresponding 
quark singlets, rifj^ and are the matrices of the Yukawa couplings. The Yukawa Lagrangian in 
Eq. (HJ) opens up the possibility that there appear tree-level FCNC. In the SM and in model I and 
model II types of the 2HDM, such FCNC at tree level are forbidden by the GIM mechanism [j|] and 
by an ad hoc discrete symmetry [^, respectively. However, tree-level FCNC are permitted in the 
general 2HDM, and this type of 2HDM is referred to as model III in the literature. 
In this model, it is possible to choose 0i and 02 in the following form 
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with the vacuum expectation values. 



<0i>=i=f M ;<02>=O. (3) 
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With this choice, the SM particles can be collected in the first doublet and the new particles in the 
second one. Further, we take Hi, H2 as the mass eigenstates h^, respectively. Note that, at tree 
level, there is no mixing among CP even neutral Higgs bosons, namely the SM one, and beyond. 

The part which produces FCNC at tree level is 

I^Y^FC = ^f.^QiLhUjR + e5Q,X02^,7i + ikM2Em + h.c. . (4) 

In Eq.(Q), the couplings ^^'^ for the flavor-changing charged interactions are 

i^h = ^neutral VcKM , 

^ch = VcKM ^neutral , (5) 

where (.n^utrai is defined by the expression 

(T> = . (6) 

and ^.neutral dcnotcd as Cn^- Here the charged couplings are the linear combinations of neutral 
couplings multiplied by Vckm matrix elements (see [ 1T0| ] for details). 

After this brief summary about the general 2HDM, now we would like to present the calculation 
of the matrix element for the B r^r^'j decay . For a general investigation of the B r^r^^ 
decay, we start with the LO QCD corrected effective Hamiltonian which induces the corresponding 
quark level process b ^ s t^t^ , given by [jTTj] 

—ACi^^ 10 ^ 

^= ^7|^W.|EC.(/^)0.(/i) + ECQ,(/i)Q.(/i)|, (7) 

where Oi are current-current [i = 1,2), penguin (i = 1,..,6), magnetic penguin (i = 7,8) and 
semileptonic (z = 9, 10) operators . The additional operators Qi, (i = 1, .., 10) are due to the NHB 
exchange diagrams, which give considerable contributions in the case that the lepton pair is r+r^ 
[|TT]]. Ci{fi) and CQ.(/i) are Wilson coefficients renormalized at the scale yU. All these operators and 
the Wilson coefficients, together with their initial values calculated at yU = rriw in the SM and also 
the additional coefficients coming from the new Higgs scalars are presented in Appendix A. 

The short distance contributions for B t^t^^ decay come from the box, Z and photon penguin 
diagrams, which are obtained from the diagrams of Fig. ([I]) by attaching an additional photon line 
either to the initial quark lines that contribute to the SD part of the amplitude, or to the final lepton 
lines, which give the so-called IB part of the amplitude. Following this framework, the general form 
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of the gauge invariant amplitude corresponding to Fig.([T]) can be written as the sum of the SD and IB 
parts 



M{B ^ T+T--f) = MsD + MlB: 



where 
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In Eqs. (^) and (|T^, and k^ are the four vector polarization and four momentum of the photon, 
respectively, q is the momentum transfer and pb is the momentum of the B meson. The form factors 
g, f, gi, /land are defined as follows [@, [T^]: 

(7(A;) |s7^(l T 75)&| B{pb)) = A^{e,,^.e*'' q^k" g{q') ± i [e*^{kq) - (e*g)P] /(g^)} , 



{^\sia^yky{l =F l5)b\B) 
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et,apaelkpq„gi{p^) ^ i e*^{kq) - {tk)q^, /i(/) \ , (13) 



(0|s7^75&|5) = - ifBPB^, ■ 
In obtaining the expressions and ([T^), we have also used 

(0|sa^,(l+75)6|5) = 
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and conservation of the vector current. Note that in contrast to AisD part of the amplitude, its M.ib 
part receives contributions from NHB exchange diagrams, which are represented by the factors F 
and Fi in Eq. (|T0|). During the calculations of these NHB contributions in model III, we encountered 
logarithmic divergences and used the on-shell renormalization scheme to overcome them. (For details, 
see [OD). 



We now examine the probability of the process B t^t^j as a function of the four momenta of 
the particles. In the center of mass (CM) frame of the dileptons r+r^, where we take z = cos 9 and 9 
is the angle between the momentum of the fi-meson and that of r^, double differential decay width 
is found to be 
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where x = 2E^/mB is the dimensionless photon energy and v = ^Jl — with r = ml/m?^. After 
some calculation, we get for the different parts of the squared matrix elements in Eq. ( [T5| ) : 
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There is a singularity in at the lower limit of the photon energy due to the soft photon 

emission from charged lepton line, while ITW^dI^ and Re{AisDM.*iB) terms are free from this singu- 
larity. It has been shown that when processes B t'^t^'j and B t^t~ are considered together, 
the singular terms in |A^/b|^ exactly cancel the 0{a) virtual correction in _B ^ r+r^ amplitude. 
But instead of this approach we prefer the one used in ref.lQ which amounts to impose a cut on the 
photon energy, i.e., we require > 50 MeV. This restriction means that we only consider the hard 
photons in the process B t^t^'^. Therefore, the Dalitz boundary for the dimensionless photon 
energy is taken as 
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with 5 = 0.01. 

Using Eqs. ([T?|)-([rSl), we get the following result for the double differential decay width 
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Integrating the angle variable, we find the photon energy distribution given by 
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We also give the forward-backward asymmetry, Aps, in B ^ t^t 7. Using the definition of differ- 
ential ApB 
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we find 
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and D{x) is given by Eq. 

Finally, we would like to discuss the r~ lepton polarization effects for the process B — > T^r~'y. 
The longitudinal polarization asymmetry of the r lepton is defined as 

{dr{SL)/dx) - {dr{-SL)/dx) 



Pl(x) 



(26) 



{dT{SL)/dx) + {dT{-SL)/dx) 
where Sl is the orthogonal unit vector for the polarization of the r lepton to the longitudinal direction 
(L) and in the CM frame of the r+r^ system, it is defined as 

f\pi\ ErPl \ 



(27) 



Here, pi and Er are the three momentum and energy of the lepton in the CM frame, respectively. 
Calculation of Pl leads to the following result 
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In order to investigate the dependence of the r lepton polarization on the model III parame- 
ters, we eliminate the other parameter, namely x, by performing the x-integrations over the allowed 
kinematical region (Eq.(|T9|)) so as to obtain the averaged lepton polarization. For the longitudinal 
component the averaged lepton polarization is defined as 



^(l-4m2/m|) ^ 
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Js ax 
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■Is d^ax 

For the process B r^r^'y, the lepton polarization has, in addition to the longitudinal compo- 
nent Pl, transverse and normal components. Since these two orthogonal components are proportional 
to the tau lepton mass, they are expected to be significant for the r+r^ channel. We shall discuss their 
effects in a more detailed paper. 

3 Numerical analysis and discussion 

To calculate the decay width, first of all, we need the explicit forms of the form factors g, f, gi and 
/i. In refs. [Q] and [jT^], they are calculated in the framework of light-cone QCD sum rules and their 
dependences, to a very good accuracy, can be represented in the following dipole forms. 
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where 

^(0) = IGeF, /(O) = 0.8 GeF, ^i(O) = 3.74 Ge\/2^ /i(0) = 0.68 G'e\/^ 
rag = 5.6 GeV , nif = 6.5 GeV , rrig-^ = GAGeV , rrif-^ = 5.5 GeV . 

In addition to these form factors, the other input parameters which we have used in our numerical 
calculations are given in table I. 
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For the free parameters of the model III, namely, the masses of charged and neutral Higgs bosons, 
mji^i, myio, m/jo, mji^^o and the Yukawa couplings (^|^'^), we use the restrictions coming from B — > 



Xg 7 decay, whose BR is given by CLEO measurement [ |T4| ] as 

BR{B^X,-i) = (3.15 ± 0.35 ± 0.32) X 10"^ (30) 

and B^ — B^ mixing [[T^, p parameter [|15]1 and neutron electric-dipole moment [|T^, that yields 
^N,ib ~ ^iid ^^ jj ~ 0, where the indices i,j denote d and s quarks, and « ^Ntt- Therefore, 
we take into account only the Yukawa couplings of b and t quarks, and ^^^^ and also Cw.rr- 
Further, in our numerical calculations we adopted the restriction, 0.257 < {Cj^^l < 0.439 due to the 
CLEO measurement, Eq.(|30|), (see [ |T0| ] for details) and the redefinition 



CU,D _ l^^F ^u,D 

Before we present our results, a small note about the calculations of the long distance (LD) effects 
is in place. We take into account five possible resonances for the LD effects coming from the reaction 
h ^ sipi ^ sT^r^, where ^ = 1, 5 and divide the integration region into two parts: 5 < x < 
1 - ((m^2 + 0.02)/mB)^ and 1 - ((m^j - 0.02)/mB)^ < x < 1 - {2mr/mBY, where m^^ = 3.686 
GeV is the mass of the second resonance. (See Appendix for the details of LD contributions). 

In this section, we first study the dimensionless photon energy dependence of the differential 
branching ratio (dBR/dx) and the model III parameters dependence of the BR and also the forward- 
backward asymmetry, ApB- The results of our calculations are presented through the graphs in 
Fig.(|t)-(|ni)- In Fig.d), we present dBR{B — ^ T^T 7) /dx as a function of x = 2E^/mB for = 
AOrrLb and C.ntt = 10 m^-, in case of the ratio \rtb\ = < 1, including the long distance contri- 

butions. Here, the differential BR lies in the region bounded by dashed (solid) curves for C^" < 
(Cj^^ > 0). We see from this figure that there is an enhancement for the dBR{B T^T^^)/dx in 
model III compared to the SM result for the Cj^^ > case, while for C^"^-^ < 0, model III predictions 
almost coincide with the SM one (small dashed curve). Fig. (|3|) is devoted the same dependence of 
differential BR, but for ^^^b = 0.1 m?, and ^^^t = "^r, in case of the ratio rtb > 1. We see that 
model III predictions for the Cj^'^ > and Cj^^ < almost coincide and they are one order larger 
compared to both \rtb\ < 1 case and the SM one. 

Fig (Q) and (|5|) show i'^^bb/^b dependence of BR for ^^^t- = , in case of the ratio \rtb\ < 1 
and r^fc > 1, respectively. In Fig.(^, BR is restricted in the region between dashed lines (solid curves) 
for Cj^^ < {Cj^^ > 0), while the small dashed straight line shows the SM contribution. In Fig. (^, 
there is a single curve since the contributions for both Cj^-^ > and Cj^^ < fit onto each other. We 
see that BR is quite sensitive to the parameter i'^bb/^b, for both \rtb\ < 1 and rtb > 1, however the 



behavior is opposite for these two cases; for \rti,\ < 1, BR is decreasing with the increasing values 
of iNbbl^b, while for rtb > 1, it is increasing. Further, BR is 2-3 orders larger compared to the SM 
result for r^f, > 1 case. For \rtb\ < 1, the enhancement with respect to the SM prediction is relatively 
moderate; nearly (30 — 40)% for C^^-^ > 0, but for C^-^^ < they almost coincide with the SM one. 

The dependence of the BR on the Yukawa coupling is presented in Fig. @ ((0)), for \rtb\ < 1 
( rtb > 1) case with ^^^^ = 40 (^^^^ = 0.1 mb). It is seen that BR is increasing with the increasing 
values of ^^^^ and this is the contribution due to NHB effects. From Fig. @, we see that the BR 
lies in the region bounded by solid lines for Cy^^ > and it is sensitive to the NHB effects, while for 
Cj^^ < 0, it is almost the same as the SM result (dashed straight line). Note that the SM prediction 
for the BR is 7.89 x 10^^ and in model III without NHB effects, when \rtb\ < 1 it is in between 
(7.83 - 7.97) X 10-9 for Cj^-^ > and (8.73 - 8.50) x lO^^ for Cj^^ < 0. When rtb > 1, upper 
and lower limits of the BR without NHB effects are (1.03 - 1.04) x 10"^ for both Cj^^ > and 
C^^^ < 0. Thus, contribution from NHB effects is seen to reach the values that are two orders of 
magnitude larger than the overall contributions for both \rtb\ < 1 and r^f, > 1, even for small values 

In Fig. (|]), the differential Apsix) is shown for ^^^^ = 40 and ^n,tt = 10 "^r, in case of the 
ratio \rtb\ < 1. Here, Apb{x) is restricted in the region between solid curves for Cy^^ > 0. It IS seen 
that the value of |ylj7B(x)| stands less than the SM one. The dashed curves represent Cy-^-^ < case 
and they almost coincide with the SM prediction for Apb{x). Fig.@ is the same as Fig. (||), but for 
rtb > 1 with ^^bb = 0.1 mfe . For this case, the sign of Apb{x) is opposite to the SM prediction and 
I^FB(a^) I is one order of magnitude smaller than the SM one. 

In Fig. ( [T0| ) we plot the Aps as a function of the Yukawa coupling for ^^^^^ = 40 mf, and 
\rtb\ < 1. ApB lies in the region bounded by dashed (solid) lines for Cj^^ < ( Cj-^-^ > ). As seen 
from Fig. (|TU|) that Aps vanishes for the large values of for C^^-^ > 0, while Cj^-^ < 0, it does 
not vanish in the given region of ^^tt- Contributions to \Apb\ from model III stand less than the SM 
ones. Fig.([TT|) is the same as Fig. ( [T0| ) , but for rtb > 1 with j,^, = 0.1 m;,. Here, contributions for 
C^^^ > and C^-^-^ < coincide and both are restricted by the solid curves. 

We present our analysis on the longitudinal component of the r lepton polarization through the 
graphs in Figs.(P^-(p3|). The dependence of on x is presented in Fig (|T^, for \rtb\ < 1 case with 
^Ntt = 10 m^. Here, Pp is restricted in the region between dashed (solid ) curves for Cy^^ < 
(Cy^^ > 0), while the small dashed curve shows the SM contribution. From this figure, we see that 
the 2HDM contributions change Pp significantly compared to the SM case for Cj^^ > 0, especially 
for the small values of x. Fig. ( [131 ) is the same as Fig. ([121), but rtb > 1 with ^^^j, = 0.1 m;,. In this 
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figure, two solid curves restrict the possible values of Pl for both Cj-' < and > and it is 
seen that both the magnitude and the sign of Pl are changed for rtt > I. 

The dependence of < > on the NHB parameter is presented in Fig. (|14[) ((|15j)), for 
\rtb\ < I (rtb > I) case with ^^^^f, = 40 (^§^bb = 0.1 nib). We see that for \rtb\ < I, < Pl > lies in 
the region bounded by dashed lines (solid curves) for Cj-^-^ < ( Cj^-^ > 0) (Fig. (p^), while in case 
of Ttb > I, the contributions for both Cy^^ < and C^-^-^ > coincide and they are represented by the 
single solid curve in Fig.([T5|). It is obvious that < P;. > is very much sensitive to the NHB effects for 
both \rtb\ < 1 and rtb > 1 cases. We note that SM prediction for < > is —0.36 and in model III 
without NHB effects, when \rtb\ < 1, it is about -0.31 for C^^^ > and -0.36 for C^^^ < 0. When 
rtb > 1, the value of < P^ > without NHB effects is about -0.36 for both Cj^^ > and Cj^^ < 0. 
Thus, the value of < Pl > without NHB effects reaches at most the SM prediction, but NHB effects 
enhance it between (5 — 100)%, even for small values of We would like to summarize our 



Parameter 


Value 


nir 


1.78 (GeV) 


rric 


1.4 (GeV) 


nib 


4.8 (GeV) 


"em 


129 




0.04 


nit 


175 (GeV) 


my/ 


80.26 (GeV) 


mz 


91.19 (GeV) 


rriHO 


150 (GeV) 


ruho 


70 (GeV) 


rriAQ 


80 (GeV) 


niH± 


400 (GeV) 


^QCD 


0.225 (GeV) 


as{niz) 


0.117 


sin^Ow 


0.2325 



Table 1: The values of the input parameters used in the numerical calculations. 

results: 

• We observe an enhancement in the differential branching ratio and branching ratio for the exclu- 
sive process B r+r^7 in the general 2HDM compared to the SM predictions. For \rtb\ < 1 
case, this enhancement is much more detectable for Cj^^ > case compared to the C^-^'^ < 
one. For rtb > 1» we see that contributions for Cj-^-^ > and Cj-^-^ < almost coincide with 
each other, and the enhancement with respect to the SM is much more sizable. 
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• BR for B T~^r~'-f decay is at the order of magnitude 10~^ (10~^) in the SM and in model III 
without NHB effects for \rtb\ < l(rtb > 1). However, including NHB exchanges may enhance 
it almost two orders of magnitude compared to the SM prediction, even for the smaller values 

• \Afb \ is at the order of magnitude 10~^ (10~^) for \rtb\ < 1 (rtb > 1) case and smaller compared 
to the SM results, which is -0.181. 

• The 2HDM contributions change and < > greatly compared to the SM case and these 
quantities are very sensitive to the NHB effects. 

In conclusion, we can say that experimental investigation of BR, Afb and Pl may provide an essential 
test for the effects of NHB exchanges and new physics beyond the SM. 
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A The operator basis 

The operator basis in the 2HDM (model III ) for our process is [ |TT| , [T7| , [T^ 

q=u,d,s,c,b 
q=u,d,s,c,b 

O5 = (sLa^tibla) (QRlBYQRf^), 
q=u,d,s,c,b 

q=u,d,s,c,b 

g 

Os = j^s^T^^a^,{m,R + msL)bf,g''>''' , 
g 

g 

= 4r^islb%)ifr), 

Qs = 4^2(^1^ r) E 

q=u,d,s,c,h 
q=u,d,s,c,b 

= 4^2(^Lb'n) E i'ii^'' 

q=u,d,s,c,b 

^2 



q=u,d,s,c,b 

q=u,d,s,c,b 

QS = ^2(-'LCT^''b%) E ifR^^^Ql) 

q=u,d,s,c,b 

Q9 = ^2(''l^"'b'R) E ifL^,M 

q=u,d,s,c,b 

Qio = ]^2i'L^"'b''n) E i<lR^,>^Ql) 

q=u,d,s,c,b 
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where a and (5 are SIJ{?)) colour indices and T^^ and Q^^ are the field strength tensors of the elec- 
tromagnetic and strong interactions, respectively. Note that there are also flipped chirality partners of 
these operators, which can be obtained by interchanging L and R in the basis given above in model 
III. However, we do not present them here since corresponding Wilson coefficients are negligible. 

B The Initial values of the Wilson coefficients. 



The initial values of the Wilson coefficients for the relevant process in the SM are []T7[] 



ctlA^w) = 

Cf-'(mvK) = 1 



4(xt - ly 2A{xt - 1; 



3 



^"^^^ - "4(x,-l)^^^^*+ 8(x,-l)3 ' 
Ct\m^) = --^Bix,) + ^ ~!^ff^'^ Cix,) - D{x,) + |, , 

C^o^'imw) = -4^iB{xt)-C{x,)) , 
sm a^r 

C^firnw) = z = l,..,10 (32) 
and for the additional part due to charged Higgs bosons are 

ctAmw) = 0, 

(mH/) = Y'F,{yt) + XY F^iyt) , 

C^{mw) = Y^G^{yt) + XY G2{yt) , 

C^{mw) = Y'H,{yt) , 

Ci'oimw) = Y'LM , (33) 



where 



Y - (^D ,?D ^ts\ 

^ — I ^N,bb "T ^N,sbT7' j ) 

nib \ ' VtbJ 

y = — Uht + iltc'^] ■ (34) 



The NHB effects bring new operators and the corresponding Wilson coefficients read as 



^N,Tri^N,ttfrnbytiQ5{yt)zA - QiizA,yt)) 



QAi^N,u) ) - 327T^mlomS,{zA,yt)e,{yt) 
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^Q2 ((^^,tt)^) 



^Q2 i^^fit) — 



Qi{zA,yt)05{yt) 



xyxt + 2za 2 Ij^^ ^^®5 ) 



-xyxtyti- 



dAn'^m^omt \Q5ixt) <di{zA,Xt) Qi{zA,Xt) 
{x - l)xt{yt/zA - 1) - (1 + x)yt) 



2pD tD 



xjyt + xtjl - yt/zA)) - 2yt 
QeQsizA) 



XfVt + 2y(a; - l)xtyt - za{xI + ©e) ^ ytjzj^ ^ ^ zaQ% 



(66 -{x- y){xt - yt))iQ3izA) + (x - y){xt - yt)zA) 

^QA^N,bb)- g4^2^2 



('AO 



r^H^atU n2n _ 9'^{iN,tt)'^'fnbmr ( Xt{l - 2y)yt , {-I + 2 cof? 6 w) {-I + x + y)yt 



cos^ OwQiivt) 



^ zh{Qi{zh, yt)xyt + cos^ Ow (-2x^(-l + a:t)|/|/| + xxtyyf - ©8^g)) \ 

cos2 6'w^ei (2:^7, 7/4)67 J' 



(35) 



-1 + 2cos^ gyt/) i/t _ xtyt _^ xtytjxy - zh) 
^biVt) QiizH^yt) 

Qb{.yt)zH 



cos2 Ow QiiVt) 

-1 + 2cos^ dw)ytZH 
cos^ ^14^67 



2xt In 



.0i(^i?,yt). 

4 _ g^mbin^xt / _ ^ , (-1 + 2x)xt , 2a;t(-l + (2 + xt)y) 



^1 ^ 1287r2m|,om2 1^ ^ 65(0;,) + y{l - Xt) 



+ 



65 (xt) 



4cos2 6li^(-l + X + y) + xt{ x + y) ^ Xt{x{xt {y - 2zh) - ^zh) + 2zh) 



cos^ 9wQi{xt) 



Ql{zH:Xt) 



yt((-l + x)xt2;ij + cos^ Owii^x - y)zH + Xt{2y{x - 1) - 2;if(2 - 3x - y)))) 



cos2 6'iv(63(^h) + - yt)zH) 



+2 {xt In 



0i(^ij,2;t) 



+ In 



xjyt - Xt)zH - Qsizn) 
(65 (xt) +y{l-Xt)ytZH 



32n^mlomtQi{zh, 1/4)65(2/*) 



(Q,{zh, yt){2y - 1) + 65(yt)(2x - l)zh) 



rho((cU ^2^ _ ^^,rr(C^,tt)^ {^b{yt)zh{yt - l)(x + y - 1) - 6i(^ft, yt){QM + ^t) iv. 



61(^^)65(^4) 



^fc€>5(yt) 
6i(zft) . 



^^^^,rr^^,ti"^6'^t /2(-l + (2 + Xt)?/) Xt{x - l){yt - Zh) 



65(Xt) 



+ 



x{xt{y - 2zh) - 4:Zh) + 2zh 

Ql{Zh,Xt) 



e'2{zh) 
(1 + x)ytZh 



+ 2 In 



^^65(3:4) 

_ei{zh,xt)_ 



xyxtyt + Zh{{x- y){xt - yt) - 65) 



^ 69 + ytZh{{x - y){xt - yt) - 65) (2a; - 1) ^ x{ytZh + xtjzh - yt)) - 2ytZh \ 
^ft 65(66 -{x- y){xt - yt)) 62(2;^) J ' 
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9^^N,rT^N,bb ( yXtVtixX^iyt - Zh) + Q(iZh{x - 2)) 



QAir'^mf^o \ 



+ 21n 



XtVt 



+ 2 In 



XtVtZh 



where 



ei(cj,A 

e5(A 

e, 



67 
eg 



= — (— 1 + y — ?/A)u; — a;(?/A + — cjA) 

= (a;* + 2/(1 - a;i))?/tt^ - xxtiyyt + (y* - l)u;) 

= Q2{'^,xt^ yt) 

= {xt{-l + y)- y)yti^ + xxtiyyt + ^^(-1 + yt)) 

= 1 — X + XLJ 

= X + A(l — x) 

= {xt + y{l - xt))yt + xxt{l - yt) 

= {yiVt - 1) - yt)zH + x{yyt + {yt - 1)zh) 

= yt{2x\l + xt){yt - 1) + xt{y{l - yt) + yt) + x{2{l -y + yt) 

+ xtil-2y{l-yt)-3yt))) 

= -x^{-l + X + y){-yt + x{2yt - l)){yt - Zh) - xtytZh{x{l + 2x) - 2y) 

+ VtiM^"^ - 2/(1 - a;)) + (1 + x)(x - 



(36) 



and 



Xt 



yt 



rat 



mt 



Zh 



Zh 



za 



mH± mj^o mlp "^^0 



The explicit forms of the functions Fi(2)(yf), Gi(2){yt), Hi{yt) and Li{yt) in eq.(Pl|) are given as 

yt{7-5yt-8y^) , y^^yt - 2) 



Fiiyt) 
F2{yt) 
Giiyt) 
G2{yt) 
Hi{yt) 



+ 



72{yt - 1)3 + 12(y, - 1)4 
yti5yt - 3) yti-^yt + 2) 

Myt - 1)^ 



I2{yt - 1)2 
l/i(-l/,2 + 5yt + 2) 
2A{yt - 1)3 

yt{yt-^) , 



+ 



%,-l)2 2(y,-l)3 

1 — Asin^Ow xyt 



w 



1 



^i(2/0 



yt - 1 (2/i - 1)2 
477/2 -7% + 38 3i/3-6|/, + 4 
18(2/t - 1)4 
1 



108(2/* - 1)=^ 



xyt 



sin^Ow 8 



+ 



y* - 1 {yt- If 



(37) 
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Finally, the initial values of the coefficients in the model III are 

ttIt sm v\Y 

C5™^(mH.) = 0, ^ = 5,...,10. (38) 

Here, we present Cq^ and Cq^ in terms of the Feynmann parameters x and y since the integrated 
results are extremely large. Using these initial values, we can calculate the coefficients Cf^^*^(/i) 
and Cgf^^^fi) at any lower scale in the effective theory with five quarks, namely m, c, d, s, b similar 
to the SM case [|lBHg ^ 0. 

The Wilson coefficients playing the essential role in this process are C|^^^^(/i), C|^^^^(/i), 
(j2HDM^^-^^ (j2HDM^^^ (j2iiDM p^j. completeness, in the following we give their explicit 
expressions. 

where the LO QCD corrected Wilson coefficient Cj^'^^^^^ (fi) is given by 

^LO,2HDM^^^ = r/16/23c2™/(^^) + (8/3) (r/^^/^S _ ^W/23^c2HDM^^^^ 

+ C2™^(mH.)EM'^' . (39) 

i=l 

and r] = as{mw) /cusil^), hi and are the numbers which appear during the evaluation pT]]. 

Cg-^-^{fi) contains a perturbative part and a part coming from LD effects due to conversion of the 
real cc into lepton pair r+r" : 



C^^if^) = Cr\f^) + Yresonis) , (40) 



where 



+ h{z, s) (3Ci(/i) + C2(/i) + SCsifi) + C,{fi) + 3C,{^) + Ce{^l)) 

- h{l, s) {4CM + ACifi) + 3C,{fi) + CM) (41) 

- ^/i(0, s) {CM + 3C4(/i)) + ^ (3C3(/i) + CM + 3CM + CM) , 
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and 



Y 



3 iiT{Vi^ T+T )mv, 
A-V = —'^ 1^ — —■ ^ 

(3Ci(/i) + C2{^l) + 3C3(/i) + C4(/i) + 3C5(/i) + CM) 



(42) 



In eq.(PDD, the functions h{u, s) are given by 

rrib 



h{u, s) 



— In In M H — - H — x 

9 fi 9 27 9 



--{2 + x)\l - x\^/^ 
9 



In 



2 arctan / , 

VX — l 



in] , for a; = ^ < 1 
forx = ^ > 1, 



h{0,s) 



27 ^ 9 



4, 4. 

- m s H — ivr 
9 9 



(43) 



(44) 



with u = The phenomenological parameter k in eq. ( p2| ) is taken as 2.3. In Eqs. (37) and (p2|), 
the contributions of the coefficients Ci(yu), Cq{ii) are due to the operator mixing. 
Finally, the Wilson coefficients Cq^ (fi) and Cq^ (fi) are given by [ pT] ] 



C^Q,(/^) = r''/''C'Q,(m^), z = l,2. 



(45) 
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Figure 5: The same as Fig.(Q), but for > 1. 
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Figure 6: BR as a function of ^§ for ^§ 



AO nib and \rtb\ < 1. 
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Figure 7: The same as Fig. (||), but for = 0.1 rub and rtb > 1. 
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Figure 8: Differential Afb{,x) as a function of x for ^^^^^ = 10 nir, ^^^b = 40 nib and \rtb\ < 1. 
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Figure 11: The same as Fig. (|T0|), but for ^^^^ = 0.1 and rtb > 1. 
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Figure 12: Pl as a function of x for ^^^^ = 40 nib and i^^rr = 10 in case of the ratio |rtb| < 1. 
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Figure 14: < > as a function of ^n,tt ^^ ^f, = 40 m?,, in case of the ratio \rtb\ < 1. 
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Figure 15: The same as Fig. (|T4|), but for ^^^^ = 0.1 and rtb > 1. 
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